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The Mixed Spin 5 = (|, 1) XXZ Ferrimagnet at Zero Temperature 
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Linked cluster series expansions about the Ising limit are used to study ground state preperties, 



viz. ground state energy, magnetization and excitation spectra, for mixed spin S = (|, 1) quantum 
ferrimagnets on simple bipartite lattices in 1, 2, and 3-dimensions. Results are compared to second- 



order spin wave theory and, in general, excellent agreement is obtained. 



O ■ I- INTRODUCTION 

Ferrimagnets are materials where ions on different sublattices have opposing magnetic moments which do not 
\ cancel in the ordered phasei. This can arise either through unequal numbers of ions on the sublattices or the ions 
having different spin quantum numbers. There is growing interest in such systems, both from fundamental physics 
and through their technological potential. Arguably the simplest such structures are bipartite lattices (A,B) with 
Sa ^ Sb , with nearest neighbour antiferromagnetic exchange coupling. The experimental discovery of ferrimagnetism 
in bimetallic chains^ has led to many studies of mixed-spin chains^^. The rare-earth nickelates R^BaNiOs 6 , which 
can be modelled by a 2-dimensional net of coupled chains^, are an example of a mixed-spin system in 2-dimensions. 



Another example is an Fe-Ni cyanide bridged network 9 , which provides a realization of the square lattice model we 
consider in this paper. 

We consider, in this paper, the cases of a bipartite chain, sqaure lattice, and simple-cubic lattice where sublattice 
A is occupied by S = \ spins and sublattice B by S — 1 spins, with nearest neighbour antiferromagnetic coupling 
between the sublattices. This is the extreme quantum limit, where quantum fluctuations will be most significant. The 



Hamiltonian is taken as 

H = J5>?SJ + + s-Sf)} (1) 

where the parameter < x < 1 represents exchange anisotropy, the coupling J > and it can be set to be 1. 

Previous work on this system includes spin- wave theory^, and a variety of systematic numerical method a?i 4 '^ , with 

most emphasis on the mixed-spin chain. Other work has included frustrating further neighbour interactions for the 

O ■ chainii and square lattice—, but we do not consider such cases here. 
CJ • ... ... 

Our approach is to use high-order linked-cluster expansions—, where quantities are expanded about the Ising limit 

as power series in x, a technique we have employed in previous studies of quantum antiferromagnets— . This allows 

the whole region < x < 1 to be studied, as well as extrapolation to the isotropic Hamiltonian x = 1. In particular 



we compute series for the ground state energy, the magnetization on each sublattice, and the dispersion relations for 
magnon excitations. 

Since we will compare the series predictions with spin-wave theory, we previde here the relevant spin-wave results 
at lowest order, for the general case with spins Si and S2. Further details, including the second-order results, are 
given in the Appendix. The ground state energy and magnetizations are 



2N ■ 

k 



Eo/NJ = -zSiS 2 /2 + ^- J2lJ (Si + S2) 2 - 4SiS 2 a: 2 7£ - (Si + &)] (2) 
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Mi >2 = S*i, 2 - ^ Y)[ ===4== ==5=f - 1] ( 3 ) 
N ^ y /(S 1 + S 2 ) 2 - 45iS" 2 x 2 7^ 

where the sum is over N/2 k-values in the reduced Brillouin zone, z is the lattice coordination number, and 7k is 
the usual factor 7k = j Yls • The first term in ©, © is the classical result, and the second term is the lowest 
correction due to quantum fluctuations. 

The two branches of magnon energies are given by 



<4 



lzJ[ y J(S 1 +S 2 ) 2 -4S 1 S 2 x 2 rt ± (S 2 - Sr)} (4) 



The most interesting feature of these is that, at x = 1, one branch is gapped, while the other becomes gapless but 
with quadratic dispersion 

k ^ (5) 

u! k ~ const, k 2 ' 

This result is well known, and is a consequence of S\ ^ S 2 - Another feature, which we have not seen commented on 
previously, is that the various quantities arc analytic at x — 1, unlike the case of the usual antiferromagnet, where 
there are square root singularities. This make the extrapolation of our series to x = 1 much more precise. 

The outline of the paper is as follows. In Section II we give (brief) details of our method and present some of the 
raw data. In Section III we present an analysis, mainly in graphical format, for 1, 2 and 3-dimcnsions. The spin wave 
results are shown for comparison. We find that many of our results agree extremely well with second-order spin-wave 
theory, and are indistinguishable on the scale of the fingures. In Section IV we present some conclusions. 



II. DERIVATION OF SERIES EXPANSIONS 

The series are obtained, as described in ref.^3], by evaluating the "proper" or "cumulant" contribution for each of 
a set of connected clusters, and then summing these, with their respective embedding factors, to obtain the series for 
the bulk system. In this way the ground state energy per site, for example, is written in the form 

oo 

E /NJ = e n x n (6) 

n=0 

where the purely numerical coefficients e n are computed exactly to some maximum order (in this work, to order 22, 
14, 12 for d = 1, 2, 3 respectively). The coefficients are given in Table I, for the case Si = ^, S 2 — 1, which is the 
only case we consider. Table II gives the corresponding coefficients for the magnetization series. We note that, apart 
from the constant term, these series are identical and hence only Mi is given. 

Series for the excitation spectra are obtained via the linked cluster method of Gelfand (see in which we write 

c(k)=^(R) e * k - R (7) 

R 

where the i(R) are a set of "transition weights" obtained from an effective Hamiltonian. The series for two branches 
of excitation spectra has been computed upto order 22, 10, 8 for d = 1, 2, 3 respectively, which involve a list of 22, 
185928, 59804 clusters, respectively. In Table II we provide our results for these, for the square lattice. In Table I, 
we also list the series for excitation spectra at k = for all three lattices. We are happy to provide other results on 
request. 
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FIG. 2: Magnetization per site on the S = \ sublattice for the linear chain (a), square lattice (b) and simple cubic lattice (c). 
The dotted and dashed lines, in each figure are the results of first-order and second-order spin- wave theory. 



III. RESULTS 



The series can be evaluated by Pade approximants or integrated differential approximantsi^, either for a particular 
x or extrapolated to x = 1. We display our results in graphical form. 



A. Ground state energy 



Figure n shows the ground state energy per site, as a function of exchange anisotropy x, for the three lattices. 

We note from this figure the excellent agreement between our series results and, particularly, second-order spin- wave 
theory. The agreement is not perfect for the chain but for d = 2, 3 the curves are indistiguishable on the scale of 
the figures. The estimated numerical error in the series is less than the width of the curve and so these results are 
essentially numerically exact. 

Figure |21 shows the magnetization as a function of anisotropy x, again for the three lattices. Let us start from the 
d = 3 case, where spin wave theory would be expected to be most accurate. We see that this is indeed the case, and 
second order spin-wave theory is barely distinguishable from our series results. At the isotropic point x = 1 quantum 



0.2 0.4 0.6 0.8 1 
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FIG. 3: Excitation energies versus k for the s — (|,1) system on the linear chain. The dotted and dashed lines give the 
corresponding first-order and second-order spin-wave theory. 

fluctuations reduce the moment from 0.5 to 0.445, an 11% reduction (actually our numerical result is 0.44450(3)). 
For the square lattice quantum fluctuations reduce the sublattice moment by about 20% at x = 1. The series results 
are close to, but distinguishable from, second-order spin-wave theory. For the linear chain the differences are greater 
and spin- wave theory does not perform well, especially approaching x = 1. The series results are still very precise, 
and give a sublattice moment of ~ 0.29 at the isotropic point, i.e. a 40% reduce due to quantum fluctuations (the 
numerical value is 0.292487(6)). We note here that Tiani& has show rigorously that this system will have long range 
order at T = 0, even in 1-dimension. First-order spin wave theory overestimates the effect of quantum fluctuations 
near x = 1, while the second-order curve shows a dip and then a sharp rise near x = 1. This is clearly a spurious 
effect. 

B. Exitated states 

Finally we turn to the excitation spectra. These are only shown for the isotropic case, x = 1. Figure |2I shows 
the two excitation branches for the chain versus k, again together with the spin-wave results. The k 2 dispersion is 
evident in both spin wave and series curves. Both bands have an approximante cosine shape. For the lower branch 
spin wave theory is quite accurate, except near k = 5, but for the upper branch first-order spin-wave theory gives a 
large underestimate. Second-order theory is better, but still under estimates the excitation energy over the entire k 
range. 

Figure^lshows the dispersion curves for magnon excitation for the square lattice, along 3 lines in the Brillouin zone. 
In this case the series and second-order spin- wave results agree very well at all points, while the first-order spin- wave 
results are too low. Note the almost flat spectrum along the zone boundary from (tt, 0) to (tt/2,tt/2). In first and 
second order spin- wave theory this is completely flat, as the k-dependence arises from 7k = ^(cosfc^ + cosk y ). The 
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(k x ,k y ) 

FIG. 4: Magnon dispersion curves for the s = (|, 1) system on the square lattice, for special lines in the Brillouin zone. The 
dotted and dashed lines give the first-order and second-order spin-wave results. 

series result shows a slight rise from (ir,0) to (ir/2,ir/2). 

The corresponding results for the simple cubic lattice are shown in Figure again for special lines in the zone. 
Agreement between series and spin- wave results is excellent, with the second order spin- wave results virtually indis- 
tiguishable from the series. 

IV. CONCLUSIONS 

We have used linked cluster series expansions about the Ising limit to investigate ground state properties of simple 
mixed spin S — (^, 1) quantum antiferromagnets, on linear chain, square and simple cubic lattices. The series are very 
regular and easily analyzable, and we believe the results are essentially numerically exact. We have calculated the 
ground state energy, sublattice magnetizations, and magnon excitation spectra, and compared our results with first 
and second order spin wave theory. As is known from previous work on Heisenberg antiferromagnets, second-order 
spin-wave theory is surprisingly accurate, even in 1-dimension. In the present case the agreement is even better, as 
demonstrated by our results above. Presumably this reflects the fact that the ground state has a finite spin. 

Our calculated excitation spectra show a quadratic dispersion at k = 0, and we have also computed the curvature 
a, defined by = ak 2 . In Table II we summarize all our numerical results at x — 1, and present a comparison with 
other methods. 

This work demonstrates that the series method is very powerful in studying mixed spin quantum ferrimagnets, and 
encourages its use in more complex mixed-spin models, including, for example, frustrated systems. We are pursuing 
a number of calculations along these lines. 
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FIG. 5: Magnon dispersion curves for the s = (^,1) system on the simple cubic lattice, for special lines in the Brillouin zone. 
The dotted and dashed lines give the first-order and second-order spin-wave results. 



Appendix 

We provide, for completeness, a summary of first and second order spin wave theory for this system. Equivalent 
treatments have been given in Ref. Iiol and elsewhere. 

The Hamiltonian for the anisotropic system is (setting J = 1) 

H = [ fi f S ™ + x ( s l S ™ + «)] + £ *f + £ S m> (8) 

(lm) I m 

where we have divided the lattice sites into even and odd sublattices, denoted by I and m respectively, and introduced 
the magnetic fields hi and h 2 . 

We firstly introduce boson operators a/ and b m via the Dyson-Maleev transformation on the two sublattices: 

1 - sublattice : sf = Si-a}a h s+ = {2Si) 1/2 a l -{2S l y 1/2 a\a l a h = {2Si) 1/2 a\, 

(9) 

m - sublattice : S z m = b\ n b m -S 2 , S+ = (2S 2 ) 1 ' 2 b\ n - (2S 2 )- 1/2 b m b m b m , S m = (2S 2 ) 1/2 b m . 

Note that this transformation is not Hermitian. In terms of the boson operators, the Hamiltonian can be expressed 
as: 

H = -7V(SiS* 2 z- h x S x +h 2 S 2 )/2 

a\ai + (zSi + h 2 ) ^ b m b m + x\J S\S 2 ^ (aib m + a\b m ) 

l m (lm) 

- ^2 aj a i b m b ™ - ^V^i^ y^^ajaiaibm/Si + a\b\j) m b m / S 2 ) ■ (10) 

(lm) (lm) 



Then, as in Ref.HJ, we introduce the Bloch-type boson operators a&, b k by a Fourier transformation: 



1/2 



iV 



1/2 



-ik-m i 



ID 



where N is the total number of lattice sites. The quadratic part of H can be diagonalized by a Bogoliubov transfor- 
mation: 



a k = a k cosh 9 k - 01 sinh 9 k , 



-at sinh 6» fe + fc cosh 6 k 



(12) 



where tanh26» fc = xj k /D, and D = + zS 2 -hi + h 2 )/{2z^/S 1 S 2 ). 

With this, one can get the ground state energy up to second order in 1/S expansion: 



E h /N = {-zS 1 S 2 + h 1 S 1 + h 2 S 2 )/2 + z^S^DC^/2 

-f {[<?-! + D D(C? - C^)] 2 + x- 2 D 2 (l - x 2 D 2 Q ){C , l 1 - C?) 2 }, 



(13) 



where D = D(h, = h 2 = 0) = {Si + S 2 )/{2^S~ 2 ), = f £ k [(l - x 2 7 ^- 2 )™ /2 - 1]. 

Setting the external magnetic field to zero (i.e. h\ — h 2 = 0), one can derive from Ea. HlHJI the ground-state energy 
E : 



o/N = -zSiS 2 /2 + z{Si + S 2 )d/4 - ~{[C-i + D 2 (d - C_!)] 2 + x- 2 D%(l - x 2 D 2 )(C- 1 - d) 2 } (14) 



where 



C B = C*= = ^E[( 1 - aa ^ IJ o 2 ) B/a - 1 ' 



Differentiating Ea. (|13|) with respect to hi, one finds the magnetization Mi: 



Mi = 



2 



N dhi 



h, =h 2 =0 



Si 



C-i C_3 — C-l 



-[D 2 {\ - x 2 )(d - + x 2 0D 2 - l)C_i 



(15) 



(16) 



2 2a; 2 (Si + S 2 ) 

Similarity, one can get M 2 , which is different from Mi by an constant S 2 — S±, as expected. 

One can also get the two branches of magnon dispersion without the external magnetic field, again up to second 
order in a 1/5* expansion: 



u±/J = ±z(S 2 -S 1 )/2 + z{S 1 + S 2 )Q 1 /2± 



D z(S 2 - Si) 



(C-i - Ci) 



(x- 2 - l)£ 2 (Q-i - Qi)(C_i - Ci) + Q1C1D 2 + Q_iC_i(l - D 2 ) 



(17) 



where Q n = (1 - x 2 j 2 D^ 2 ) n ^ 2 

For x = 1, we recover the second-order spin- wave results of Ivano\iii. It is interesting to study the asymptotic 
behaviour for the quantities C n near x = 1. For the case Si = S 2 , the leading order correction term to C n near x = 1 
ialA const x (1 — a; 2 ) 1 / 2 , while for case ^ S2, it is wasy to show that the asymptotic behaviour for the quantities 
C n near x = 1 is: 



C n (x) = C n (l) + -[C n _ 2 (l) - C„(l)](l - x 2 ) + 



(18) 



8 

For linear chain (ID), square lattice (SQ) and simple cubic lattice (SC), C n at x = 1 for (Si, S2) = (1/2, 1) are 

ID Ci = -0.29097039335, C_i = 0.60977301072, C_ 3 = 5.3812664598 (19) 
SQ Ci = -0.13351109261, C_i = 0.21847549012, C_ 3 = 1.3753265565 (20) 
SC Ci = -0.08401123644, C_i = 0.11629596999, C_ 3 = 0.5647562506 (21) 

(22) 
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TABLE I: Series coefficients for the ground state energy per site Eq/NJ Magnetization Mi, and two branches of magnon w^/J 
at k = for the linear chain, the square lattice and the simple cubic lattice. Nonzero coefficients x r up to order r = 22 for Id, 
order 14 for square lattice, and order 12 for simple cubic lattice are listed. 



r 


Eo/NJ 


Mi 


w-(k = 0)/J 


cu + (k = 0)/J 






linear chain 









-5.00000000xl0 _1 


5.00000000xl0 _1 


1.00000000 


2.00000000 


2 


-2. 50000000x10^ 


-2. 50000000x10^ 


-2.00000000 


-1.00000000 


4 


1.04166667xl0~ 2 


-5.90277778xl0~ 2 


1.12500000 


2.31250000 


6 


2.56799769xl0~ 2 


2.01786748xl0 _1 


3.45052083 


-5.81163194 


8 


-1.40457216xl0~ 2 


-7.01659922xl0~ 2 


-3.34128599 


1.92529744xl0 +1 


10 


-6.52895221 xl0~ 3 


-1.51931908xl0 _1 


-1.24692841 xl0+ 2 


-7.43345195 xl0 +1 


12 


1.20231578xl0~ 2 


1.67134766X10" 1 


9.95763258xl0+ 2 


3.08779448 xlO+ 2 


14 


-9.68145478xl0~ 4 


6.73415017xl0~ 2 


-3.98363592 xl0+ 3 


-1.34245308xl0+ 3 


16 


-1.00328363xl0~ 2 


-2.45737502 x 10" 1 


7.07316949xl0+ 3 


6. 03441420 xlO+ 3 


18 


6.99312641 xl0~ 3 


8.82808138xl0~ 2 


6.98417467xl0+ 3 


-2. 78235446 xl0+ 4 


20 


5.40374090 xl0~ 3 


2.41826850xl0 _1 


5.10317554xl0+ 4 


1.30865902xl0+ 5 


22 


-1.04177497xl0~ 2 


-2.78535289xl0 _1 


-1.73937095 xl0+ 6 


-6.25428515xl0+ 5 






square lattice 









-1.00000000 


5.00000000X10" 1 


2.00000000 


4.00000000 


2 


-2.00000000xl0 _1 


-8.00000000xl0~ 2 


-1.90000000 


-1.40000000 


4 


-5. 22222222 xl0~ 3 


-1.25320988xl0~ 2 


-1.08660714xl0~ 2 


2.88214286xl0~ 2 


6 


-7.77432120xl0~ 4 


-2. 94020395 xl0~ 3 


-6.48931742 xl0~ 2 


-3.05217039xl0~ 2 


8 


-3.41477688xl0~ 4 


-1.37128294xl0~ 3 


-9.96489680 xl0~ 3 


-1.95162662xl0~ 2 


10 


-9.77862982xl0~ 5 


-5.45189634xl0~ 4 


-6.07211572xl0~ 3 


5.18453462xl0~ 4 


12 


-4.20558132xl0~ 5 


-2.71359073xl0~ 4 






14 


-1. 75132661 xl0~ 5 


-1.33427314xl0~ 4 










simple cubic lattice 









-1.50000000 


5.00000000X10" 1 


3.00000000 


6.00000000 


2 


-1.87500000xl0 _1 


-4.68750000xl0~ 2 


-2.46428571 


-2.03571429 


4 


-3.89229911 xl0~ 3 


-4.79023836xl0~ 3 


-2.21734002 xlO -1 


-2.12222866xl0 _1 


6 


-1.24575187xl0~ 3 


-1.91047285xl0~ 3 


-1.50472245 xl0 _1 


-1.23920805 xl0 _1 


8 


-4.21416977xl0~ 4 


-9.04898794xl0~ 4 


-5.75189462 xl0~ 2 


-5.89627745 xl0~ 2 


10 


-1.66838367xl0~ 4 


-4.59923079xl0~ 4 






12 


-7.78978452xl0~ 5 


-2.59993308xl0~ 4 
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TABLE II: 

COs(nk x ) COS 


Series coefficients for two branches of magnon dispersion u> (k x ,k y ) = 
(mk y )]/2 on the square lattice. Nonzero coefficients a,k. n ,m up to order k 


= 10 are 


n v k 

listed. 


cos(mk x ) cos(nk y ) + 


(k,n,m) 




(k,n,m) cifc,n,m 


(k,n,m) a k , n 


m 


(k,n,m) 





LO' /J 



( 0j 


0. 


0) 


2.000000000 




< 8 


2) 


3.713762427x 10" 


-3 


f 6 


3, 


3) 


-8 349006559 x 10" 


-3 


( 8 3 5) 


-2 354734575x10" 


-3 


( 2, 


o, 


0) 


-4.000000000x10" 


i 


(10, 0, 


2) 


1.371857188x10 


-3 


( 8 


3, 


3) 


-2.461673445x10" 


-3 


(10, 3, 5) 


-1.492768067x10" 


-3 


( 4, 


0. 


0) 


6.325595238x10" 


2 


( 4 2 


2) 


-5.812500000x10" 


-2 


(10 


3, 


3) 


-5.509512597x10" 


-4 


( 8 2 6) 


-1.177367287x10" 


-3 


( 6, 


0. 


0) 


-5.757990909x10" 


3 


( 6 2 


2) 


-8.531165228x10" 


-3 


( 6 


2, 


4) 


-1.252350984x10" 


-2 


(10 2 6) 


-1.124360790x10" 


-3 


( 8, 


0. 


0) 


-1.507778153x10" 


4 


< 8 2 


2) 


1.296551584x 10" 


-3 


f 8 


2, 


4) 


-4.921900666x 10" 


-3 


f 8 1 71 


-3 363906535 x 10" 


-4 


(10, 


0. 


0) 


4.027910703x10" 


4 


flO 2 


2) 


-7.669154691x10" 


-5 


(10 


2, 


4) 


-1.312852923x10" 


-3 


(10 1 7) 


-5.883082637x10" 


-4 


( 2, 


1, 


1) 


-1.000000000 




( 4, 1, 


3) 


-7.750000000x10 


-2 


( 6, 


1, 


5) 


-5.009403935x10" 


-3 


( 8, 0, 8) 


-2.102441585x10" 


-5 


( 4 , 


1, 


1) 


7.619047619x10" 


2 


( 6 1 


3) 


-2.053764395x10" 


-2 


( 8 


1, 


5) 


-3.914975780x10" 


-3 


(10, 0, 8) 


-9.061290707x10" 


-5 


( 6, 


1, 


1) 


8.413191563x10" 


3 


( 8 1 


3) 


1.132685760x10" 


-3 


(10 


1, 


5) 


-1.582142621x10" 


-3 


(10, 5, 5) 


-3.071949277x10" 


-4 


( 8, 


1, 


1) 


3.139563846x10" 


3 


(10, 1, 


3) 


-6.338716677x10 


-4 


( 6, 


0, 


6) 


-4.174503279x10" 


-4 


(10, 4, 6) 


-5.119915461x10" 


-4 


(10, 


1, 


1) 


2.939474184x10" 


3 


( 4, 0, 


4) 


-9.687500000x10 


-3 


( 8, 


0, 


6) 


-8.955696061x10" 


-4 


(10, 3, 7) 


-2.925665978x10" 


-4 


( 2, 


0. 


2) 


-5.000000000x10" 


1 


( 6, 0, 
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TABLE III: Comparison of some numerical estimates obtained by different methods for the ground-state energy Eo/NJ, the 
magnetization Mi, and the energy gap uj + (k = 0), and a at x = 1. 



method 


Eo/NJ 


Mi 


cj + (k = 0) 


a 






linear chain 






series (present work) 


-0.7270467(10) 


0.292487(6) 


1.7591(6) 


1.66(3) 


DMRG 3 


-0.72704 


0.29248 






QMC£ 






1.75914(1) 


1.48(4) 


1st order SWT 


-0.718228 


0.195113 


1 


2 


2nd order SWT 


-0.730420 


0.316344 


1.67556 


1.5218 


3rd order SWTJA 


-0.727161 


0.293884 










square lattice 






series (present work) 


-1.2065125(10) 


0.40206(1) 


2.5775(8) 


1.80(4) 


1st order SWT 


-1.20027 


0.390762 


2 


2 


2nd order SWT 


-1.20731 


0.401293 


2.52798 


1.87255 






simple cubic lattice 






series (present work) 


-1.693375(7) 


0.44450(3) 


3.505(10) 


1.69(8) 


1st order SWT 


-1.68903 


0.441852 


3 


2 


2st order SWT 


-1.69371 


0.444025 


3.45069 


1.95157 



